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We consider the combined effect of a gap and a Zeeman interaction on the helical Dirac fermions 
which exist on the surface of a topological insulator. Magneto-optical properties, the magnetization, 
Hall effect and the density of states are considered with emphasis on the particle-hole asymmetry 
which arises when a subdominant Schrodinger piece is included along with the dominant Dirac 
part of the Hamiltonian. When appropriate, we compare our results with those of a single valley 
gapped graphene system for which Zeeman splitting behaves differently. We provide a derivation of 
the phase offset in the magnetic oscillations brought about by the combined effect of the gap and 
Schrodinger term without requiring the semiclassical Onsager quantization condition. Our results 
agree with previous discussions based on semiclassical arguments. 
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I. INTRODUCTION 

A topological insulator (TI)ii^ is a material which is 
insulating in the bulk but posses a metallic spectrum of 
surface states. These are expected to have high mobil¬ 
ity as they are protected by the topology of the band 
structure. The surface states of Bi2Se3 have a single 
Dirac cone^ centred at the T-point of the two-dimensional 
( 2 D) surface Brillouin zone. Other materials can have a 
different odd number of such points. As an example, 
samarium hexaboride has three^. Such a new phase of 
matter comes about through a band inversion caused by 
a strong spin-orbit coupling. At present, these systems 
are extensively studied due to their novel physics and 
their possible applications in devices such as quantum 
computing^ and photonics^. An important observation 
is the phenomenon of spin-momentum locking^ estab¬ 
lished by spin- and angular resolved photoemission spec¬ 
troscopy (ARPES). This shows that the in-plane compo¬ 
nent of the spin is perpendicular to the in-plane momen¬ 
tum. A gap can also be introduced at the Dirac point 
through doping of the TI surface with magnetic (time- 
reversal-breaking) impurities. In the work of Chen et 
al^, a gap of A ~ 7 meV was observed through ARPES 
measurements on a sample of (Bio,ggMno,01)2803 with 1 % 
Mn substitution on the Bi site. Another avenue to cre¬ 
ate a gap is to make the TI ultra thin with a distance 
between the top and bottom surface of the order of the 
extent in space (perpendicular to their plane) of the sur¬ 
face states^^— . Tunnelling between the top and bottom 
surfaces gaps the electronic spectrum. Of course, in this 
case, both surfaces can contribute to a particular phe¬ 
nomenon and a sum over both cones (which have gaps of 
opposite sign) is needed as has been discussed recently by 
Yoshimi et alf^i and Zhang et aljiS in the context of the 
quantum Hall effect. Neupane et alf^^ find that varying 
the quantum-tunnelling gap in ultra-thin films leads to a 
modulation of the in-plane spin texture. While the spin- 


momentum locking remains, the in-plane component of 
spin can itself be decreased with decreasing thickness. 
Other such studies include work by Tahir et ali^ on the 
oscillations expected in the quantum capacitance of thin 
films. For simplicity, we will consider a single gapped 
Dirac cone; however, the results applicable to thin films 
can be obtained by assembling two such cones with the 
appropriate relationship between their respective gaps. 

Particle-hole asymmetry is an important feature of the 
surface-fermion band structure of a TI. This arises from a 
subdominant quadratic-in-momentum Schrodinger term 
with mass m in addition to the dominant linear-in- 
momentum Dirac term with Fermi velocity vp. The rel¬ 
ativistic piece has its origin in the strong spin-orbit cou¬ 
pling and involves the Pauli spin-matrices sigma^ and 
ay. This is distinct from the case of graphene where the 
matrices exist in pseudospin space and the energy lev¬ 
els are degenerate in spii>i^. With the application of a 
magnetic field perpendicular to the plane, the Landau 
levels (LLs) are split by the Zeeman field and form two 
sublevels with a definite Sz component of spin 
This contrasts with the case of the TI for which only the 
N = 0 level is fully spin polarized in while the av¬ 
erage Sz for all the other levels is small in comparison. 
For the presently studied TI, the Dirac magnetic energy 
scale at one Tesla is an order of magnitude larger than the 
Schrodinger scale. However, at B = 36 T (for example) 
the two scales differ by less than a factor of two. Even 
when B is small, the Schrodinger term can lead to novel 
effects such as the splitting of the optical absorption lines 
seen in the real part of the AC conductivityii. Thus, this 
term must be included in any complete theoretical treat¬ 
ment and is important when considering experimental 
data. In particular, optical data gives valuable informa¬ 
tion on the dynamics of the surface charge carriers in a 
TI as demonstrated in the experimental work by LaForge 
et alji^ and Orlita et alj^^ among others. 

Our paper is organized as follows: in Sec. H, we specify 
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our low-energy Hamiltonian involving Schrodinger-Dirac 
kinetic energy and the gap term which does not 
commute with the rest of the Hamiltonian. The eigenen- 
ergies and eigenfunctions which arise in the presence of 
a magnetic field are worked out; from these, the aver¬ 
age value of Sz for each level is determined. We also 
consider a related system written in pseudospin space. 
In this case, the wave-functions are simultaneous eigen¬ 
states of the Hamiltonian and Sz. Thus, each LL is 
fully spin polarized in the z direction. In Sec. HI, we 
begin our discussion by considering the integrated den¬ 
sity of states up to energy Wmax (as measured from the 
Dirac point). We emphasize the effect of particle-hole 
asymmetry in a TI and contrast this with the case of a 
single valley of a graphene-like system which shows no 
asymmetry in the absence of a gap. We also discuss 
the modifications brought about by a finite A and by 
Zeeman splitting. Section IV contains formulas for the 
magneto-optical conductivity. Both the longitudinal and 
the transverse Hall conductivities are considered. Again, 
particle-hole asymmetry is emphasized. We find a com¬ 
parison with the single-valley graphene-like system to be 
helpful. Circular dichroism is discussed. In Sec. V, we 
examine the magnetization of the metallic surface states. 
Through the Streda relation22i2i, we relate its derivative 
with respect to chemical potential (/r) to the quantiza¬ 
tion of the Hall plateaus. As expected, we find complete 
agreement between this and our conductivity formula. 
We find that the quantization remains half-integral (i.e. 
±1/2, ±3/2, ±5/2,...) which is characteristic of relativis¬ 
tic fermions (with due consideration for the degeneracy 
factors). The value of /r at which a new step occurs is 
changed with the value of Schrddinger mass, A and Zee- 
man splitting. In contrast to graphene, Zeeman interac¬ 
tions do not introduce additional stepsi^ between those 
present when the effect is neglected. It is important to 
realize that when the Dirac term in our Hamiltonian is 
dominant, this describes a TI and the Hall quantization is 
relativistic. When the Schrddinger piece dominates, the 
Hall quantization is non-relativistic and the Hamiltonian 
applies to the usual spintronic materials. In Sec. VI, we 
consider the magnetic oscillations and in particular, em¬ 
phasize the phase offset in the usual semiclassical expres¬ 
sion. This quantity is 1/2 for non-relativistic Schrddinger 
particles and 0 for relativistic Dirac fermions. In the pres¬ 
ence of a finite mass and magnetic gap, we find that the 
offset is changed to 7 = — A(1 ± g)/{2mvp), where g 
is the Zeeman strength. This reduces to the relativistic 
result for either A = 0 or m —>■ 00. The phase offset 
depends on the Zeeman term through the coupling coef¬ 
ficient g. Except for this factor, our quantum mechani¬ 
cal results agree with a previous result^ which relied on 
semiclassical arguments2^“— . Fuchs et al/2i start with 
the Bohr-Sommerfeld quantization condition suggested 
by Onsager— and include a correction in the band ener¬ 
gies due to the magnetization M (fc) induced by the exter¬ 
nal magnetic field B of the form —A4{k) ■ B. They apply 
this to a two band model of gapped graphene and explic¬ 


itly show that the phase offset (7) in the magnetic oscil¬ 
lations is zero in this case even though the Berry phase 
around the cyclotron orbit is not tt. Because of the gap, it 
is instead equal to ttWc [1 — A/where Wc is the wind¬ 
ing number and is a topological invariant referred to by 
Fuchs et al2i as the topological part of the Berry phase. 
They establish that 7 = 1 /2—IVc/2 = 0 in this case and it 
is only the topological part of the Berry phase which en¬ 
ters this quantity. Thus, the —Wc/2 exactly cancels the 
Maslov index contribution of 1/221. Following this line 
of reasoning, Wright and McKenzie^ extend this semi¬ 
classical method to the case when a Schrddinger mass is 
introduced as well as a gap. In this case, the result for the 
phase offset is given by their Eqn. (26). To lowest order 
in the Schrddinger mass, this reduces to A/{2mvp). We 
derive this result without resorting to any semiclassical 
arguments, and further find its generalization to include 
the Zeeman term. 


II. SURFACE STATE HAMILTONIAN 

In the presence of magnetic dopants, the helical sur¬ 
face states of a TI are given by the Bychkov-Rashba 
Hamiltonianllill 

H =-^+ hvFika:^y - kya^) + A&z, ( 1 ) 

where ax, ay and dz are the usual Pauli-spin matrices 
and fik is the momentum measured relative to the F- 
point of the surface Brillouin zone. The first term is the 
Schrddinger piece for describing electrons with effective 
mass m. The linear-in-momentum term describes mass¬ 
less Dirac fermions which move with a Fermi velocity vp- 
The last term accounts for the magnetic impurities which 
open a gap of 2A in the band structure. For Bi 2 Te 3 , band 
structure calculationsHiH predict: vp = 4.3 x lO^m/s, 
and m = O.OQrUe where rUe is the bare mass of an 
electronll. Solving Eqn. m yields the energy dispersion 

E± (fc) = ^ ± ^{hvpky ± AA (2) 

2m 

A schematic plot of the surface-state band structure is 
given in Fig. [TJ Figure UKa) shows the gapless spectrum 
which we contrast to frame (b) where a finite A is in¬ 
cluded. It is important to note that we are mapping 
a continuum Hamiltonian onto a finite Brillouin zone. 
Equation [T] permits the valence band to bend back to¬ 
ward the zero energy axis; this is unphysical for TIs. As 
a result, an appropriate momentum cutoff must be ap¬ 
plied. 

To examine the magnetic response of such 
systemsiSi^i"— , we must consider the effect of a fi¬ 
nite magnetic field B oriented perpendicular to the 
surface (z). In the Landau gauge, the magnetic vector 
potential, defined through B — V x A, may be written 
as A = {0,Bx,0). We account for the magnetic field by 
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FIG. 1. (Color online) (a) Gapless spectrum of topological 
surface states, (b) The effect of including magnetic dopants 
on the low-energy band structure. In both cases, significant 
particle-hole asymmetry is observed as the mass term causes 
the conduction band (blue) to narrow while the valence band 
(red) fans out. 


applying the usual Peierls substitution to the momentum 
{hk hk + eA) in Eqn. (IT|). Our Hamiltonian becomes 

H = — [{-iOxf + {-idy + eBx/Kf] A Aaz 
+ hvF[{—idx)d'y — {—idy + eBx/K)ax\- (3) 


A Zeeman interaction can be included by adding 
— {l/2)gsiiBBaz, where gs is the coupling strength (Ri 8 
for Bi 2 Se 3 ^) and p-s = eh/{2me) ~ 5.78 x 10“^ meV/T 
is the Bohr magneton. Next, we define the raising and 
lowering operators 


and 


= 


72 L 


-ikx + 


a; -I- So 


/2 

‘‘B 


a = 


IB 


72 I 


ikx 


s -I- So 


(4) 

(5) 


where Ib = \/Ti/{eB) and sq = kyPg. On a Fock state 
(|A^)) of the harmonic oscillator Hamiltonian, these op¬ 
erators have the property: 

a|Af) = 7fV|fV-l), (6) 


and 


(7) 


With these definitions, the Hamiltonian reduces to 


/ 


H = 


ml^ L 


a^a - 


Z 


-Tlvf-j—C l 


—avp-—a 
Ib 


ml'^ 


a7 - 


-Z 


( 8 ) 

where Z = A — gsgBB/2. Equation ([8]) yields the LL 
dispersion 


Sn,s = 


EqN + s\ 2NEf + 


A-^{l+g) 


1 2 


iV = l,2,3,. 


^(1+5)-A 


N = 0 
(9) 


where Ei = hvp \/eB/h, Eq = HeB/m, g = ggm/{2m,e) 
is the renormalized Zeeman coupling coefficient, and s = 
± for the conduction and valence band, respectively. It 
is important to note that in a sum over s, the N = 0 level 
will only contribute once. Unless otherwise noted, we will 
take Ei/VB = 10.4 meV/7T, Eq/B = l.I meV/T and 
H = IT. We note that Eqn. ([5]) is identical to that given 
previously by Shen et al4^i^. However, these authors 
are only interested in the limit of spintronic materials 
for which the Schrodinger term in Eqn. m is dominant. 
Here, we treat the opposite limit where the Dirac term 
is dominant (requiring a band cutoff). Thus, there is 
no overlap between the two studies and the works are 
complimentary. 

Since our Hamiltonian is written in the spin basis, the 
eigenstates are comprised of spin-up and -down ampli¬ 
tudes. Indeed, the wave-functions of Eqn. ([8]) are 


where 


\Ns) 




( 10 ) 


F — 


-S\ 


'1 , A-{l+g)Eo/2 


2{£n,+ - EoN) 

0 


fV=l,2,3,... 
iV = 0 


( 11 ) 


is the spin-up component and 


C 


I _ 

N,s — 



A - (1 + g)Eo/2 
® 2{£n,+ - EoN) 

1 


iV = 1,2,3,... 

N = 0 

( 12 ) 


gives the spin-down amplitude. We immediately see 
that the N = 0 LL is entirely populated by spin-down 
electrons^ (i.e. ^ = 0). To see the z component of 

spin of the remaining levels, we compute the average 

(7) = (A^l ^7 |A^) • 


7 |7V) = TFTT |fv + 1). 


(13) 
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This gives 

; f ^ ^[A-(l + g)i^o/2] 

(«-)=2) ^2ElN+[^-{l + g)E^/2f 

y -1 7 v = o 

(14) 


where ■C = i for the K and K' points of the hexagonal 
Brillouin zone. Applying the same techniques as before, 
the LL energies are^l^— 

_ / —^gsUBBcr + S\/ + 2NEf, iV = 1, 2, 3,... 

-^A-^g,f,BBa, N = 0 

(16) 


Clearly, an spin polarization is quickly lost for N ^ Q. 

Examining the N = 0 level [Eqn. ([9])] reveals that for 
A = 0, the LL sits at positive energy (i.e. fo,+ = (1 + 
(/)£’o/2). Eor finite A, the level remains above zero as 
long as {l+g)Eo/2 > A. For A > {l+g)Eo/2, the zeroth 
level is situated at negative energy. Also, in contrast 
to graphene (where spin is a good quantum number), a 
finite Zeeman term does not split the levels but simply 
renormalizes their energy. As opposed to graphene, the 
N ^ 0 levels are not simply shifted by Egs^BB/2. The 
zeroth level is translated in the usual way since it has a 
definite spin-down polarization. A schematic plot of the 
low-energy LLs in a TI is shown in Fig. [2] Particular 

A=0 (l+g)Eo/2>A (l+g)Eo/2<A 



(a) 


(b) 


(c) 


where cr = ± for spin up and down, respectively. The 
corresponding wave-functions are 


and 


\Ns)j, = 


—iAN,s |A^ — 1) \ 

Bn,s \N) ) 


\Ns)j,, 


—iAN,s I A') \ 


respectively, where^ — 


and 



A^O, 

A = 0, 



A^O, 
A = 0. 


(17) 


(18) 


(19) 


( 20 ) 


Unlike the Tl, spin is a good quantum number and thus 
a Zeeman interaction simply shifts the spin-up (-down) 
levels down (up) by a constant amount. 


FIG. 2. (Color online) Schematic illustration of the Landau 
level energies for a Tl with (a) A = 0, (b) (1 -|- g)Eo/2 > A 
and (c) A > (1 -I- g)Eo/2. As A increases, the zeroth level 
begins at positive energy and decreases (becoming negative 
when A > {l + g)Eo/2). 


attention is given to the relative location of the A = 0 
level. The various LL energies are given by the coloured 
circles (blue for conduction levels, red for valence levels, 
and purple for A = 0) and are overlaid on the B = 0 
band structure. 

Throughout this paper, we will be interested in com¬ 
paring the results of a Tl with the familiar case of gapped 
graphene; in particular, the graphene result at the K- 
point. In the presence of an asymmetry gap A, the 
graphene Hamiltonian written about a single valley ^ in 
the sublattice basis becomes 


H, = 


hv{^kx + iky) 


hv{^kx - iky) 

-A 


(15) 


III. DENSITY OF STATES 

We begin our analysis with a comparison of the inte¬ 
grated density of states. In a magnetic field, the density 
of states is given by a sum of 5-functions which peak at 
the various LL energies. For a TI, 


A(.) = f 


S{uj - £o,+) + ^ 5 (a; - £n,s) 


iV=l 

s=± 


( 21 ) 


where £i\f,s is given by Eqn. ®. Similarly, for graphene 


= X 


OO 

iV=l 
s=± 


( 22 ) 
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where is given by Eqn. (HU). To compare the two 
results, we define the total integrated density of states 
up to energy w^ax as 


r^ma.x 

d(Wmax)= / N{uj)dLU. (23) 

Jo 

A plot of /(wniax) as function of the cutoff energy Wmax 
is shown in Fig. [S] For graphene, the steps are of dou¬ 
ble weight due to spin-degeneracy. Figure |3Ka) shows the 



FIG. 3. (Color online) Integrated density of states for 
graphene at the A-point and a TI (a) without and (b) with a 
gap of A = 0.5 meV. (a) For A = 0, particle-hole symmetry 
is present in graphene. For a TI the valence and conduc¬ 
tion levels are asymmetric, (b) For finite A, the symmetry 
in graphene is removed. A Zeeman interaction spin-splits the 
graphene levels and simply renormalized those of a TI. 

A = 0 results for graphene at the A-point (solid black) 
and a TI (dashed red) in the absence of Zeeman split¬ 
ting. For graphene, the N = 0 LL is pinned at zero 
energy and particle-hole symmetry is present for all N. 


For the TI, the N = 0 level sits at positive energy and 
clear particle-hole asymmetry exists for the N ^ 0 levels. 
The valence levels sit closer to Wmax = 0 than the corre¬ 
sponding conduction levels. In Fig. [3Kb), a finite gap is 
included and the effect of Zeeman splitting is explored. 
For graphene with gg = 0 (solid black), the gap moves 
the N = 0 level to negative energy at the A-point and 
breaks the particle-hole symmetry. With a finite Zeeman 
term (dash-double-dotted green), the steps split into two 
(emphasized by the green circles). This is characteristic 
of systems with spin-degeneracy. For a TI (dashed red), 
the gap shifts the energy levels and a finite Zeeman in¬ 
teraction (dash-dotted blue) does not split the steps but 
further renormalizes their energy. A plot of the LLs for 
A = 0 and varying gs is shown in Fig. For graphene 

(a) Graphene K 
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FIG. 4. (Color online) A = 0 LLs in (a) graphene and (b) 
a TI with and without Zeeman splitting, (a) In graphene, a 
Zeeman interaction spin-splits all the LLs and shifts the spin- 
down levels (red) up while the spin-up LLs (blue) decrease in 
energy, (b) For a TI, Zeeman effects do not split the levels. 

at A [frame (a)], the gs = 0 levels are symmetric. A 
finite Zeeman splitting causes the spin-down levels (red) 
to increase in energy by gsgBB/‘J while the spin-up levels 
(blue) decrease by the same amount. This splitting does 
not depend on N and creates a series of spin-polarized 
levels. The TI [frame (b)] is particle-hole asymmetric 
and a Zeeman term causes the conduction levels (pink) 
to increase in energy while the valence levels (purple) 
decrease. Unlike graphene, the A 0 levels are not 
spin-polarized and the splitting is not uniform; for larger 
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TV, the renormalization becomes less prominent. 


IV. MAGNETO-OPTICAL CONDUCTIVITY 


In the one-loop approximation, the optical conductiv¬ 
ity aapiS^) is given by the familiar Kubo formula 




ih 


/m.s' — fN,s 


2ttI^ £]\[ s — £ 

B N,M=0 


M,i 


s,s'=± 


Ns\j^\Ms')(^Ms'\j 0 \Ns 


+ Sm,s' — 


^N,s 


ih/ (2r) ’ 


(24) 


where fn,s is the Fermi function for state n in band s, 
r is the relaxation time and ja = eva = {e/K){dH/dka) 
is the current operator. In the zero temperature limit, 
fn can be replaced by the Heaviside step function Q{pL — 
£n,s)- The necessary velocity operators are 


.. k . ^ h — a 

Vx = —kx + ypUy = i - — - 1 - VpOy, 

m Tn y/21 b 


(25) 


and 


V 


V ~ 


h 

m 



- VfCTx 


h + a 
m y/ 2 lB 


VF^x- (26) 


Calculating the appropriate matrix elements, we obtain 
the real and imaginary parts of the longitudinal conduc¬ 
tivity 



r 

(o -|- £ m , s ' ~ £ n , s )'^ + 

X [J^{Ns; Ms')Sn^m-i + J^{Ms';Ns)Sm,n-i] 


and 



T £m^s' £n,s 

{ft + £m,s' — £n,s)'^ + 

X [T{Ns; Ms')Sn,m-i + Ns)Sm,n-i] 

respectively, where F = 7i/(2r) and 


T{Ns-Ms') = 


X 


ni 

^M,s'^N,s 


Eq 



'^N,s 


+ Vn + 


(29) 

n 2 


Likewise, the real and imaginary parts of the transverse 
Hall conductivity are 



fl -|- £m,s' ~ £n,s 
"" {n + £M,s'-£N,s)^ + r^ 

X [J^(A^s; Ms')5n,m-i — E{Ms'; Ns)5m,n-i], 


Im 


^xy (^) 

e'^/h 



E 

N,M^0 

s,s^=± 


fM,s' — In,: 
£n,s — £. 


M,s' 


-F 


(31) 


(H + £m,s' - £n,s? + r 2 

X [F{Ns\Ms')5n,m-i — E{Ms'\ As)(5M.Ar-i]- 


For the gapped graphene results, the reader is referred 
to Eqns. (12)-(15) in Refi^. Note that an overall minus 
sign should be included in Eqns. (14) and (15) of Refi^. 
To obtain the appropriate limit, Ago = 0 and A 2/2 = A. 
As previously mentioned, a momentum cutoff must be 
applied to the band structure to ensure that the valence 
band does not bend back across the zero energy axis. 
When a magnetic field is applied, a corresponding cutoff 
on N must be applied to ensure that £n,- < 0 . 

We begin our discussion by consider the effect of finite 
A on the absorptive part of the longitudinal response. 
This is shown in Fig. [S] for a TI. Recra;a;(H) is plotted for 
/r = 0, g = 0 and three values of A: 0 (solid black), 0.5 
meV (dashed blue) and 1 meV (dash-dotted red). The 
inset shows a schematic plot of the lowest LLs for the 
different A’s. The lowest optical transition is marked by 
the arrows. For A = 0, the V = 0 LL sits at Eq/2 (for 
g = 0). Therefore, the lowest transition occurs between 
the TV = 1, s = — and N=0 levels. For finite A less 
than Eo/2, the N = 0 level is still at positive energy but 
its magnitude has decreased; this causes the lowest ab- 
’sorption peak to move down in H. As A becomes larger 
than Eo/2, the N = 0 level moves to negative energy and 
the first transition is now between the zeroth level and 
the N = 1 LL of the conduction band. Now, the fre¬ 
quency of the first transition continues to increase with 
A as £o^f is pushed further down in energy. Higher op¬ 
tical transitions are present and occur in pairs. This is 
a signature of the particle-hole asymmetry as the energy 
for the £n,- to £n±i,+ transition is not the same as the 
to £n,+ transition. For graphene, the two sets of 
’split higher-energy peaks in Fig. [5] would each coalesce 
into a single line. It is the presence of the Schrbdinger 
magnetic energy Eq which sets the scale for the split¬ 
ting. This energy increases linearly with magnetic field 
and inversely with decreasing m as discussed by Li and 
Carbotteii. 

Next, we explore the results for finite chemical poten¬ 
tial and 5 = 0. In Fig. [BJ the effect of positive and neg- 
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FIG. 5. (Color online) n — 0 longitudinal conductivity in a 
TI for varying A. The energy of the first absorption process 
is schematically illustrated by the arrows in the inset. As 
A increases, the zeroth LL decreases in energy. This shift is 
evident in the onset frequency of the first optical transition. 


ative fi is examined. In frame (a), the results of gapless 
graphene at the if-point are shown; due to the particle- 
hole symmetry, the response is identical for ±/i. The op¬ 
tical transitions which lead to this set of absorption lines 
are shown in Fig. [7Ka) where the arrows are color-coded 
to correspond to Fig. Ha). Note that the positions of the 
En.- LLs are at the negative of the £n,+ levels. For posi¬ 
tive /i, the five arrows on the left (black) apply; while, for 
negative it is the five green arrows on the right which 
are relevant. In both cases, each transition has a one- 
to-one correspondence with the other set. Figure Hb), 
displays the results for a TI when A = 0. Here the finite 
Schrodinger term breaks particle-hole symmetry. This is 
clear from the transitions shown in Fig. Hb) . In this case, 
even for A = 0, the negative energy levels do not mirror 
the positive energy set. The N = 0 level is no longer at 
zero energy but rather has been pushed to positive en¬ 
ergy Eo/2. The £n,+ energy is also larger than 
The black arrow between iV = 0 and iV = 1+ (which ap¬ 
plies to the lowest line for positive fi) is longer by Eo/2 
than the arrows between the N = 1~ and N = 0 levels 
which is the first transition when fi is negative. Thus, 
in Fig. Hb), the first peak of the negative fj, response 
(purple) is lower in H than the corresponding positive /r 
peak (black). The results for a gapped TI are shown in 
frame (c) of Fig. [5] with the corresponding optical transi¬ 
tions shown in Fig. He). Again, an obvious asymmetry 
is present between the ±/x regimes. The asymmetry is 
now much larger than that shown in frame (b) for two 
reasons: the finite gap adds asymmetry; the larger value 



FIG. 6. (Color online) Longitudinal conductivity in (a) gap¬ 
less graphene, and (b) a gapless and (c) gapped TI for positive 
and negative jj.. (a) For graphene, particle-hole symmetry en¬ 
sures the ±/i results are identical. For a TI [(b)-(c)], the 
results are different which emphasizes the asymmetry of the 
LLs. The shading under the dashed blue curve is for empha¬ 
sis. Note the missing peak in this case. 


of chemical potential enhances it is well. An additional 

























































effect which needs to be emphasized is that one of the ab¬ 
sorption lines for /i < 0 (blue) is missing in the first set 
of split peaks shown at higher energy. Of the four peaks, 
the first, third and fourth apply to all three curves (/J. = 0 
and p, = ±14 meV), but the second is only present in the 
black and red curves. The transitions which give rise to 
these plots are shown in Fig. [7] where the missing tran¬ 
sition is indicated by a blue x. It cannot occur because 
the N = 1~ level is unoccupied. 
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B=1T, gs=0 
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FIG. 7. (Color online) Absorption processes which give rise 
to the conductivity curves of Fig.[6l The arrows are coloured 
to correspond with the curves in Fig. |6] and the location of p 
is given by the dotted lines. For p = —14 meV in the gapped 
TI [blue arrows in frame (c)], the transition from Ti,_ to £ 2 ,+ 
is forbidden [note the missing absorption line in Fig. [^c)]. 

The absorptive part of the transverse Hall conductivity 
is also of interest and is shown in Fig. |S] for the same pa¬ 
rameters as Fig. ETc). When the optical transition is from 
the level in the valence band to the (TV — I)**' LL, the 
response is negative; while the transition from the 
to {N ± 1)*^ level is positive. This will have important 
ramifications on the circular dichroism. The shading un¬ 
der the dashed blue curve for negative ^ again helps to 
emphasize the missing peak around ^ 36.6 meV which 
occurs for ^ = 0 (black) and ^ = 14 meV (red). The 
other three peaks exist for all /i considered here. Note 
that, expect for the lowest energy peak in the dashed 
curve (blue) and in the dash-dotted curve (red), all other 
peaks would not exist in the pure Dirac limit. Their 
existence depends on the presence of the subdominant 



FIG. 8. (Color online) The absorptive part of the trans¬ 
verse Hall conductivity for a TI with parameters set to match 
Fig. Etc). For transitions from N to N — 1, the response is 
negative; for N to N + 1, the conductivity is positive [see 

Fig-EDc)]- 


Schrodinger contribution in our Hamiltonian, even if it 
is small. This represents a qualitative difference between 
the physics of the two cases. 

In the context of TI thin films, both surfaces become 
important; the two surfaces are marked by an opposite 
sign of A. In Fig. El the longitudinal conductivity for 
A = ±5 meV is shown. The combined result that would 
be measured in an optics experiment is given by the solid 
black curve. Frame (a) corresponds to charge neutrality 
(/i = 0) while finite chemical potential is considered in 
frame (b). Aside from an energy shift in the peak associ¬ 
ated with N = 0, we note that the spectral weight for the 
A —>■ A ± 1 transitions (doublets of higher energy peaks) 
is markedly different depending on the sign of the gap. In 
the dash-dotted blue curve, the first peak of the doublet 
has the largest spectral weight while it is opposite for the 
dashed red curve. 


A. Circular Dichroism 

The response to circularly polarized light is given by 
o'xxi^) iicFxyiS^) for right- and left-handed polarization, 
respectively. Therefore, the absorptive part is deter¬ 
mined by 


Re(T±(0) = Recra;a;(D) lui(Jxy{p)- (32) 

This is readily evaluated by utilizing Eqns. and m- 
A plot of the circular dichroism is given in Fig. [TU] for 
a gapped TI. Frames (a) and (b) show the response 



















































































9 



FIG. 9. (Color online) The absorptive part of the longitudinal 
conductivity for a TI with A = ±5 meV for (a) /i = 0 and 
(h) fj. = 7 meV. 


to right- and left-handed polarized light, respectively. 
Right-handed light selects out the transitions between 
the level in the valence band and the {N — 1)**' LL 
in the conduction band. Left-handed light selects the 
TV to -I- 1 transitions. Right circularly polarized light 
shows a low-energy absorption peat at ^ = 0 (solid black) 
which moves to lower energy in the dashed blue curve for 
/i = —14 meV [frame (a)]. Such a peak is missing for 
/i = -1-14 meV (dash-dotted red). Instead, this peak is 
present for left circular polarization [frame (b)] while the 
other two values of /i have no such absorption feature. 
The pairs of peaks shown at higher energy in both frames 
are present for all three values of for right-handed light. 
For left-handed polarization, the higher peak in the pair 
is present for all /r considered here, the lower one is miss- 



FIG. 10. (Color online) Conductivity response to (a) right- 
and (b) left-handed circularly-polarized light, (a) For right- 
handed light, only transitions from N to N — I are present, 
(b) In the response to left-handed polarization, only the N to 
N + 1 processes occur. 

ing for /r = 14 meV (dashed blue). This can be traced to 
the forbidden transition shown in Fig. [7l[c) (marked by 
the blue x). 

V. MAGNETIZATION OF THE METALLIC 
SURFACE STATES 

To discuss the magnetization of the surface states, we 
turn to the grand thermodynamic potential 

n{T,n) = -T f iV(w)ln (l-h dw, (33) 
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where T is the temperature (in units of fcs), and N{oj) 
is given by Eqn. (EU. The quantities discussed herein, 
depend on a i? derivative of the grand potential; for con¬ 
venience, we add (/r/2) J^^N{uj)duj to Eqn. (1^ . The 
integral of the density of states over all energies gives the 
total number of states (which must be independent oi B). 
Thus, this term will not contribute to the magnetization 
{—did/dB). At zero temperature, Eqn. (1551) becomes 


= [ (u! — N{uj)du! + f {u! — fi)N{uj)du! 

J-OO ^ 2 / Jq 

+ ^ [ N{uj)duj. (34) 
2 Jo 

For a gapped TI, 

/*0 pff 

/ N{uj)duj= / N{ijj)duj-\-—T, (35) 

Jo J-oo ^ 


where 


-1, (l+g)Eo/2<A 
T=<( 0, (1 + 5)Eo/2 = A . 

1) (1 + 9)Eo/‘2 > A 


(36) 


Therefore, 


fl(/x) = f {u! — fi)N{uj)dui-\—[ ujN{uj)duj. 
Jo 2/l J-oo 

(37) 


The final term does not depend on /r and gives the vac¬ 
uum contribution. This will simply provide a constant 
background to the /r dependence of the magnetization. 
Keeping only the /i dependent pieces, the first two terms 
of Eqn. (|57|) give 


n{p) = 


eB 


sgn(£o.-s) i£o.+ - m) 0 (ImI - 1^0,-si) 0 (sgn(£:o,H 


+ “ M) 0 “ ^N, + ) 


iV=l 


- ^ {£n- - /i) 0 {£n- - 9) 


N=1 


(38) 


where again, we require all the s = — states to be neg¬ 
ative. Note that 0(0) = 1/2 as only half the ^-function 
situated at w = 0 is integrated. 

The slope of the magnetization is of particular inter¬ 
est as it is related to the quantized Hall conductivity 
through the Streda relation^S [an = edM/d^). To see 
the quantization of the slope, note that 


dfi 


sgn(fo,-s)0 (ImI - l^o.-sl) 0 {sg^{£o,+)9) 

T 

N^l 

(39) 


The slope of the magnetization is then 


dM 

dfi 


d dQ e 
dB d^i h 


sgn(£:o,+)0 {\^i\ - |£:o,-s|) 0 isgn{So,+)9) 


- V + E ® - ^^.+) - ® 


N^l 


(40) 


A plot of the magnetization calculated from Eqn. (I58ll 
through M = —dfl/dB is shown in Fisf. ITlT al. The cor¬ 
responding slope [see Eqn. dini)] is given in Fig. db). 
Here, g is taken to be 1 and the effect of varying A is 
emphasized. A saw-tooth oscillation pattern is present in 
M(/i) with the location of the teeth sitting at the various 
LL energies. The derivative of M(/r) gives the quantiza¬ 
tion of the Hall conductivity. A half-integer quantiza¬ 
tion is present which is characteristic of Dirac systems. 
The Hall conductivity is an = /h)v with filling fac¬ 
tors V = ±1/2, ±3/2, ±5/2,.... For zero gap (solid black 
curve), the Hall conductivity at /i = 0 has v = —1/2. 
This results from the finite value of the N = 0 level^. 
As A is increased, the location of the N = 0 step moves 
lower in g,. At /r = (1 ± g)Eoj2 — A, the step occurs at 
zero chemical potential; for larger A, the small |/r| value 
of an is given by ^ = 1/2. Again, asymmetry is seen 
between the negative and positive g regimes. These re¬ 
sults have been verified by taking the DC limit (H —^ 0) 
ofRecra;y(H) [see Eqn. (150)) ]. 

To compare this with gapped graphene, the reader is 
referred to Refi^. The quantization of the Hall conduc¬ 
tivity is given by their Eqn. (15). To obtain the gapped 
graphene result, one must take Ago = 0 and A 2/2 —>• A. 
In Fig. [T51 the slope of the magnetization is shown for 
graphene with A = 2 meV, Ei/'/B = 25.64 meV/VT 
characteristic value for graphene) and B — IT. The 
upper two frames show the contribution from the K- and 
AT'-point, respectively. The lower frame shows the total 
result. Note that the spin-degeneracy is included. For a 
single spin species, the Hall response is characterized by 
a half integer filling factor. The location of the N = 0 
LL is symmetric between the two valleys (sitting at —A 
for K and A for K'). As a result, the total system near 
/r = 0 is insulating (aji — 0) since the individual val¬ 
leys contribute equal and opposite edge channels. Here, 
a Zeeman interaction (dashed purple curve in the upper 
frame) splits the steps in the conductivity into two. 


VI. MAGNETIC OSCILLATIONS 


We now wish to explore the quantum oscillations which 
exist for low B fields. This can be done by extracting the 
oscillating part of the density of states, using that to find 
the grand potential and taking the appropriate deriva¬ 
tive (M = —dVl/dB). To leading order in 1 /(tou|.), the 
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FIG. 11. (Color online) (a) /i dependence of the magnetization 
for a TI with varying A. (b) The corresponding slope of 
the magnetization (which is related to the Hall conductivity 
through the Streda relation: dM/dfj.= {l/e)aH)- In all cases, 
a half-integer quantization is present; however, the step to 
-1-1/2 can be tuned from positive to negative by increasing 
A. 


magnetization is (see Appendix) 


mLM « - 


e — A^) r, 


2TTkh 


1 + 


2^mvi 


sin (27rfca:i), 
(41) 


where 


Xi 


1 - 


- A^] 


2El 





Id (meV) 


FIG. 12. (Color online) The slope of the magnetization for 
gapped graphene at (top) the A-point, (middle) the A'-point 
and (lower) the total response. The step which results from 
the A = 0 level is at /r = =pA for the A- and A'-point, 
respectively. Thus, the total system is insulating for —A < 
jj, < A. Zeeman (dashed purple curve in the top frame) splits 
the steps into two. 


Comparing Eqn. to the customarjj^^ 


Xi 


hA{ii) 

2TTeB 


(43) 


the coefficient of the (x 1/B term in Eqn. (l42l) is in 
fact related to the area of the cyclotron orbit: 


A{fj,) = nkp 





(44) 


The remainder is the phase shift which is independent of 
B. It is. 


A(l + g) 

2mVp 


(45) 


In the pure Dirac limit (m — >■ oo, A —0), A(/i) 
reduces to the correct valu o^^i^^ of /{h^Vp). For 
finite m and A = 0, we obtain the correction of 
— [TTfi'^/{?i‘^Vp)][^/{mvp)] to A{^)^. Clearly, the gap¬ 
less limit has a phase shift of 0 associated with a Berry’s 
phase of tt— T he amplitude of the quantum oscilla¬ 
tions [in units of —e/(27rfch)] is 


Am — 


{jA - A^) 



An 


2^mvp 


(46) 
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which properly reduces to the results of Sharapov et 
al^ when to —>• oo, i.e. Am = (see their 

Eqn. (8.10) in the pure limit), and to the results of Tabert 
and Carbotte^ when A = 0, i.e. Am = ^[1+ ^/(2 tou|.)] 
[see their Eqn. (49)]. Wright and McKenzie^ applied 
the semiclassical quantization method of Onsager— . This 
was augmented by including a first correction in the band 
structure for the energy shift due to the magnetic re¬ 
sponse of the system (as discussed by Fuchs et al2i). 
They obtained a phase offset in the magnetic oscillations 
which, to within a sign that can be traced to the chirality 
of the Hamiltonian, reduces to our result [Eqn. (HSl) ] for 
large to and when the Zeeman interaction is neglected 
(i.e. g = 0). 


VII. CONCLUSIONS 

The simplest description of the surface states of a 
TI is a single Dirac cone centred at the T-point of the 
2D surface Brillouin zone. The presence of a subdom¬ 
inant quadratic-in-momentum Schrodinger term to the 
purely relativistic linear-in-momentum Dirac Hamilto¬ 
nian is also an important feature; it reshapes the perfect 
graphene-like cone into an hourglass shape with the nar¬ 
rowing of the conduction cone and outward fanning of 
the lower cone (valence band). Particle-hole symmetry is 
lost which can have important ramifications on physical 
properties such as the density of states, magnetization 
and the optical properties (both AC and DC). 

Doping the surface of a TI with magnetic atoms breaks 
time-reversal symmetry and creates a gap of 2A at the 
Dirac point. Alternatively, a slab can be made thin 
enough that the top and bottom surface states hybridize; 
consequently, both become gapped. The sign of the gap 
is opposite on the two surfaces. We have studied the 
magneto-optical response as a function of energy of such 
systems including as well the effect of a Zeeman inter¬ 
action. Particular emphasis is given to the particle-hole 
asymmetry brought about by the Schrodinger mass term 
associated with the non-relativistic piece of the Hamilto¬ 
nian. In a finite magnetic field, both the Zeeman interac¬ 
tion and the gap can modify this asymmetry. Comparing 
the integrated density of states /(wmax) to Wmax in a TI 
to a single valley of graphene highlights the important 
differences between the two cases. Namely, particle-hole 
asymmetry and the very different effect of Zeeman split¬ 
ting. For a TI, Zeeman coupling does not split the steps 
of /(wmax) into two spin-polarized substeps displaced by 
a constant amount as in graphene. Instead, it shifts the 
onset of the various steps by a level-dependent amount. 
For a single valley of graphene, the gap does produce 
asymmetry; but, this will not be seen if the two valleys of 
opposite chirality are superimposed. Since the Hamilto¬ 
nian for a TI involves real spin (as opposed to pseudospin 
in graphene and related materials such as M0S2 and sil- 
icene), only the N — 0 level is fully spin-polarized^^. All 
other LLs are found to have a much reduced value of the 


average Sz ^ (?i/2)[A — (I -I- g)Eo/2]/{Ei'/2N). This is 
to be contrasted with the pseudospin case where all lev¬ 
els are fully Sz polarized with Sz = Ehj2 depending on 
whether the LL has moved up or down in energy due to 
Zeeman. 

The dynamic longitudinal magneto-optical conductiv¬ 
ity for a TI also displays asymmetry between positive 
and negative values of the chemical potential fi. Even 
for g — 0, the absorption line associated with the £ 1 ^- 
to fo,+ transition can be eliminated and replaced by a 
fo,-(- to £i^+ line by manipulating the gap value. There 
are peaks in the imaginary part of the transverse Hall 
AC conductivity which do not exist for the pure Dirac 
case. These features do not depend on the subdominant 
non-relativistic term being large. They correspond to the 
interband optical transitions and come as pairs of peaks 
(one positive the other negative). They translate into 
new absorption peaks for circularly polarized light. 

While the starting formulas for the magnetization and 
the DC limit for the real part of the Hall conductivity 
are quite different, the values of the Hall plateaus that 
ultimately result are identical and agree with those for 
the pure relativistic case. Although, the values of chem¬ 
ical potential at which a new step appears are different. 
They depend on the presence of a non-relativistic piece in 
the Hamiltonian. For example, the transition form —1/2 
to 1/2 no longer corresponds to zero chemical potential. 
This transition is shifted to the position of the iV = 0 
LL which is given by (1 -I- g)Ao/2 — A. We emphasize 
that this quantity remains finite when both the gap (A) 
and Zeeman term {g) are zero. The shift of Eq/2 is di¬ 
rectly due to the non-relativistic Schrodinger term in our 
Hamiltonian. Even if this term is small, it provides a 
qualitative modification of the physics as compared with 
the pure relativistic case. Note that including a Zeeman 
term and the gap further shifts the energy at which the 
Hall conductivity transitions from —1/2 to 1/2. It also 
depends on the sign of the gap. For A > 0, the shift 
is to lower energies while for A < 0, it is toward higher 
energies. 

We consider the quantum oscillations that arise in the 
magnetization. Unlike previous work, our considerations 
do not involve any semiclassical arguments based on On- 
sager’s quantization condition for cyclotron orbits. Here 
we proceed from the grand potential and use a Poisson 
formula to obtain the low-field limit {B —>■ 0). We find 
that when the gap is finite and a Schrodinger mass term 
is also included, there is an offset in the phase shift asso¬ 
ciated with the magnetic oscillations. To lowest order, it 
is given by 7 = — A(1 -|- g)/{2mvp) which reduces to zero 
when A = 0 or TO —>■ 00. In both these cases, 7 = 0 as ex¬ 
pected for Dirac fermions. We also see a dependence on 
the Zeeman splitting (g). In this formula, vf is the Dirac 
Fermi velocity. Except for a sign due to the chirality of 
the Hamiltonian, this agrees with Eqn. (26) in Refj^^ 
when only the leading order in 1 /to is retained. These 
authors employed a semiclassical approximation to ob¬ 
tain their results. We note that this phase offset remains 
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even though the quantization of the Hall plateaus is rel¬ 
ativistic (±1/2, ±3/2, ±5/2,...). A new expression for 
the amplitude of the magnetic oscillations is given which 
properly reduces to that of gapped graphene when the 
subdominant Schrodinger term is dropped and to that 
previously found in Ref?^ when the Schrodinger mass is 
included but the gap is set to zero. 


with i = 1,2. This generalizes Eqn. (8) of Ref.^ to 
include a gap and Zeeman splitting. For Eq and 

m —>■ 00, 


Xi 




UJ 

mvp 


bEo 

2mvp 


ujbEo 


(A.6) 


where 
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Appendix 


and,thus 


[A-(l±g)Ao/2]^ 

El 


(A.7) 


bEo 

2mvp 


2mVpEo 


A2-EoA(1±5) 


Eq 

T 


(l±5)^ 


(A.8) 


To obtain the magnetic oscillations, return to Eqn. (ED 
and express it as 


N{uj) = 


eB d 
h doj 


0 (w — £o,+) ± ^ ^ 0 (w — £n,s) 


N=1 

s=± 


(A.l) 


The oscillating part of the density of states can be ex¬ 
tracted by applying the Poisson formula 

1 l‘00 

E ("^-2) 

N=1 

OO pco 

±2E / E{x)cos{2'iTkx)dx. 

Following arguments given by Suprunenko et al^, this 
gives, 

± [ 0 (a; ± |fo,-|-|) ± 0 (w — |fo,-|-|) ~ Q{u} — CCniin)] 


Xi 


E ^sin(27rA:a:i) 

^' TTK 


k=l 


± 0(w — Wmin) 


X2 


—sin(27rfca;2) 

^' TTK 


k=l 


where 


__El_EofEo \_A 

2Eo 2El[2^^^^ ^ 


(A.3) 

(A.4) 


and 

= TT + TT ± (-)* 





2Elu: , , 

fl + g 


El > 

K 2 

Eo) 


(A.5) 


The last term in the above expression goes like Eq (x B 
and will drop out in the limit B —^ 0^. This leaves 


bEo 

2mvp 


A2 ^ A(1±5) 


2Ef 


2mvp 


(A.9) 


where the last term is constant in B and will therefore 
contribute a constant phase to the quantum oscillations. 
We now have 


Xi 


1 


1 - 


r 2 A 21 ^(1 ± 9) 


1 - 


To lowest order in l/(mvp) 


Xi 


1 


1 - 


(jj 


[02^ - A^] ± 


A(l±g) 

2mvp 


(A.IO) 


(A.ll) 


As we do not allow the valence band to bend back toward 
the zero energy axis, and we focus on w > 0, X 2 can be 
ignored and the oscillating part of the density of states 
is entirely determined by 


Aosc(w) = |go,+ |) E ^sin(27rfca:i). 

^ fc=i 

(A.12) 


All the information about the quantum oscillations is 
contained in the first term of Eqn. (ED. The oscillat¬ 
ing part of the grand thermodynamic potential is thus, 

f^osc(/^) “ f g)A qsc. (A.13) 

Jo 


where 


Nosciu})- ^ 


(A.14) 
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and 

= 0(0. - Ifo.+ I) 

Therefore, 


sin(27rfcxi) 

irk 


(A.15) 


Consolidating these results, 
e I 


Ksc(m) 


2TTkh 


'0 


eB 

^oscitA = J 0(w - |-Eo.+ |)sin(27rfcxi)dw. 

° (A.16) 

The magnetization is given by M = —dQ./dB. Using 


Trky I „ , 
—tvtCOS I zirk 
Ef 


sin |^27rfc 

2Ef 

1 


2Ef 


where 


dxi 


dB 2ElB 
the magnetization can be written as 


1 ^ 

y - A^] , 

(A.17) 

a = 


mVp 




mVp 


yjy + IS? mv% 

- A^] , 


dy. (A.24) 


(A.25) 


TTkh L 


\A-{l+g)Eo/2\ 


sin(27rA:a;i) (A.18) 
duj. 


1 ^ 

^ 2 

[w^ — A^] cos(27rA:a;i) 

mVp 



Trk 

" W 

In the limit of interest {Eq —>■ 0), the lower bound of 
integration can be replaced by A. To solve the integral, 
define 


y = 

1 ^ 

[w^ - A- 


mvp_ 



r 3w^ 

A2 1 

dy = 

2w-, 

+ 9 


L mvp mvp j 

J 

dy 


cllo — ■ 

3uj‘^ 

A2 ■ 

2CJ- TT 

+ 2 


mVp 

mVp 


dto 


(A.19) 


For m —>■ oo, y « so, to a first order correction 

in 1/m, 


1 - 


\/2/ +A2 


mv. 


[UJ^ - A2] 


Expanding further, we obtain 


UJ 


\ 


1 + 


Vy + ^ 


and 


Therefore, 




1 + 


mvp 


■A2, 


+ A^ 


(A.20) 


(A.21) 


(A.22) 


2uj — 


3u}^ 


mv^p mVp 


A" 2^/y + A2 


1 + 


\/2/ +A2 


(A.23) 


and 


d = 


A(l+g) 

2mVp 


(A.26) 


Next, define y = aa; to obtain 
r-l 


MLid) 


2'Kkh 


I 


I y'ax + A2 mvj 


X {sin (aa; + d) — axcos [ax + d)} , (A.27) 


where a = irka/El, 6 = 2Trkd, and again we work to 
lowest order in l/(mvp). To proceed, consider 


h= (sin [ax + d) — axcos [ax + d)} 

Jo 

. r r\ dx 

= / sin(aa; + o)— , = 

Jo i/a; + A^/a 

X , d 


dx 


^/x + A2/c 


, ,_ida;—sin (aa; + d) . (A.28) 

Jo v/^TAV^ dx ^ ^ ’ 

Integrating the second term by parts and defining w = 
x + IS? j a and b = 6 — alS? we obtain 


^l+A^/a 

Ii= sin (aw + b) 

Ja^/o 


3 ^ A'^ja 


l2i/uJ 2w^/2 

sin [a + 6) 


du! 




a 


(A.29) 


To solve the remaining integral, we use the definition of 
the Fresnel sine and cosine integrals: 

S{z) = J sin dt, (A.30) 


and 


C{z) = J cos dt, (A.31) 


respectively. We make use of the fact that 
sin (aw + b) = sin(aw)cos6 + cos(aw)sin5. 


(A.32) 
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We have 


^i+aV« ^ 

/ sin [au! + o)—j= 

lA^/a 



i/l+A2/a 

,_ cosb 

i/l+A2/a 

,_ sin 6 . 

(A.33) 


Once evaluated at the limits, we will have Fresnel func¬ 
tions with arguments proportional to y/a. We are in¬ 
terested in the limit a —?► oo. We use the asymptotic 
expansions of the Fresnel integrals for a —>■ oo: 


and 

+ (A.35) 

Therefore, Eqn. (IA.33I) is proportional to 1/a and is 
negligible in the limit of interest. We now return to 
Eqn. (IA.29I) and consider 


^l+AVa duj 

/ sm (auj + 0 ) 

lA^/a 


= -^sin(aa; -I- b) 

vw 


^3/2 

h-aVq 

A^ja 


— 2-\/27ra 


S I ]lj sin5 — C 



^l+A^ja 

yfAFJa 

(A.36) 


Applying the Fresnel expansions for a —> oo, this term is 
zero. We arrive at the simple result 


sin(a -I- (5) 
a/1 -I- A^/a 


(A.37) 


Next, we consider 

nl 


l 2 = {sin (aa: -I- <5) — axcos {ax -I- (5)} dx 
Jo 

= [ sin (ax + S)dx— f xdx—sin (ax + S) . 

Jo ' Jo dx ^ 

(A.38) 


Integrating the second term by parts, we obtain 
2 - 1 

I 2 = -cos {ax + S) — sin {a + 6) . (A.39) 

a 0 

In the limit a —>■ 00 , 

I 2 ~ —sin(a -I- d). (A.40) 

Combing Eqns. (jA.37l) and (IA.40|) . the magnetization be¬ 
comes 


Mo'sc(d) « - 


-b 


a 


2Trkh ly/a + mvp_ 


Written in terms of the original variables. 


« - 


e {fjJ — A^) 

2TTkh /r 


1 -b 


2fj,mvp 


sin(a -b S). 

(A.41) 


sin (27rfca;i), 
(A.42) 


where 


Xl 


1 - 


iL 

mur 


- A^] 


2 f;/ 


-b 


A(l + g) 

2mvp 


(A.43) 
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